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i-^ ' Abstract 

1) , We propose new equations of motion under the theory of the Brownian motion to connect 

the states of quantum, diffusion, soliton, and periodic localization. The new equations 
are nothing but the classical equations of motion with two additional terms and the one 
of them can be regarded as the the quantum potential. By choosing a parameter space, 
various important states are obtained. Further, the equations contain other interesting 
'^ ■ phenomena such as general dynamics of diffusion process, collapse of the soliton, the 

nonlinear extension of the Schrdinger equation, and the dynamics of phase transition. 
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1 Introduction 

The description of microscopic world governed by quantum mechanics is quite different from 
the one of microscopic world with classical mechanics [H |2]. Whereas, it is an important 
issue to develop a framework that includes both of classical and quantum physics for the 
purpose to study mixed dynamics [U [3l [5] . Formally, one needs to handle large number of 
quantum particles to organize macroscopic world. The method treating quantum many-body 
system is accurate but time-consuming to search a numerical solution unless one can find a 
good approximation for individual circumstances. On the other hand, extension of classical 
equations of motion by adding quantum potential which yields the Schrodinger equation is 
much simpler to involve classical-quantum dynamics [6l [71 [8]. One natural motivation of 
the quantum potential approach is to explain macroscopic quantum phenomena [9], [10] . For 
the macroscopic level, classical equations of motion may be applicable even for quantum 
phenomena. In this paper, we propose more general equations under the theory of the 
Brownian motion as the origin of quantum potential to connect some famous equations of 
classical and quantum physics. As a result of this connection, quantum, diffusion, soliton, 
and periodic localized states can be derived by solving the general equations of motion. 

Basic concept is to generalize the idea of the Nelson's stochastic equations [11]. These 
are based on the Langevin equations and the Fokker-Planck equations to describe non- 
differentiable process. They can lead the Schrodinger equation by assuming some conditions 
such as the universal relation and the specific form of the mean acceleration. We generalize 
the Langevin equations so that both of the Nelson equations and the Schrodinger equation 
are extended. Including the nonlinear Schrodinger equation (NLS), we provide several new 
equations that contain interesting solutions. 

The extension of the Nelson equations was already attempted in the literature. In the 
paper [12], Davidson generalized the equations to remedy the universal relation: original 
relation is h = 2mD and modified one is ^ = a/1 — 2smD, where m is a particle mass, D 
is the diffusion coefficient, and s is a free parameter. This modification not only reconstruct 
a hypothesis of the relation but also it is possible to make physical continuation to the 
diffusion equation by assuming 1 — 2s < 0. Moreover, as is known that the limit h ^ 
recovers classical physics, choices of the parameter as s = 1/2 — h'^/2m'^D'^, s = 1/2, and 
s > 1/2 represent quantum, classical, and diffusion equations, respectively. Meanwhile, there 
is another extension of the Nelson equations with a different motivation [13]. We found a 
way to put them together in a framework of generalized equations of motion. 

The purpose of our generalization is to seek a proof of quantum potential in a similar 
way of the work [2]. With the modified Nelson equations, new static solutions are easily 
obtained. For instance, the soliton solution in our model is one of a simple prediction because 
it can be distinguished from other soliton solutions (see e.g. [15]) by comparing the density 
function. Besides, we also predict new dynamics such as a generalized diffusion process by 
inventing new diffusion equations. In several points, the existence of quantum potential in 
the origin of the Brownian motion can be examined. 

This paper is organized as follows. In section 2, we present generalized equations of 
motion with the generalization of Nelson's equations including a brief review. In section 3, 
the Schrodinger equation is derived from the equations of motion with the extension to a 



Hilbert space. Section 4 deals with the diffusion equation by choosing a parameter space of 
the equations proposed in section 2. Section 5 leads the soliton equation to exhibit new type 
of soliton solution. In section 6, we obtain the periodic localized solution that exists in other 
parameter space of generalized equations of motion. Section 7 is devoted to the summary. 

2 Modified Nelson equations 

The Langevin equation can describe the Brownian motion with a random movement. Nor- 
mally, this equation computes a motion of differentiable process. Suppose x(t) is a stochastic 
process, mean forward derivative is given by V^x{t) = \iinAt^+o{{x(t + At) — x{t))/At) and 
mean backward derivative is V_x(t) = \iinAt^+o{{x(t) — x{t — At)) /At), where () repre- 
sents mean value. For differentiable process, mean derivatives are equal T>j^x{t) = 'D^xit). 
However, there are many important processes which are not differentiable. To extend the 
equation to apply non-differentiable processes, independent equation is set for mean forward 
and backward velocities individually. Basic procedure is formalized by Nelson in |llj . 

Let us consider kinematics of Markoff process. The position x{t) of the Brownian particle 
satisfies 

d^x = b^(t,x(t))dt + dw, d^x = b_(t,x(t))dt + dw, (1) 

where b+ (b_) is the mean forward (backward) velocity and w describes the Wiener process 
which is Gaussian with mean zero. For the Brownian motion, it is usually set (dwidwj) = 
2D6ijdt. The Langevin equations provide the space-time evolution of b±. Here we assume 
modified differential equations: 

d+b^ = si{d+ - d^)b^ - S2d+d+b^ - (3vdt + K{x{t))dt + dB{t), 

d-b+ = -si{d+ - d^)b+ + S2d-d.b+ + Pvdt + K{x{t))dt + dB{t), ^ ^ 

where m^v is the force of friction, mK is a body force {irtK = — V^), S is a Wiener process 
representing the residual random impacts, and dB is Gaussian with mean zero. With free 
parameters Si and S2, the Nelson equations are extended. 
Suppose b±(t,x) can be expanded in Taylor series, we get 

rf+6_ = dtdtb. + {{b+dt + dw) ■ V)6_ + -{dw)i{dw)jdidjb_, 

d^b+ = dtdtb+ + {{b^dt + dw) ■ V)6+ - -{dw)i{dw)jdidjb+, (^) 

d+d+b^ = {dw)i{dw)jdidjb_, d_d-b+ = {dw)i{dw)jdidjb+, 

where order [dt)"^ is neglected. Taking the average of Eqs. ([2]), our generalized Langevin 
equations read 

{dt + b+-V + DV^ - si{b+ - b_) ■ V - 2D{si - S2)V^)b^dt = -(3vdt + K{x{t))dt, 

(9t + 6„ ■ V - DV^ + si{b+ - 6_) ■ V + 2D{si - S2)V'^)b+dt = (3vdt + K{x{t))dt. ^ ' 

Dividing by dt and summing the two equations, we obtain 

dtV + v-Vv -u-Vu- DV^u + 2siu ■ Vu + 2D{si - S2)V^w = K, (5) 



where v = (b_|_ + fo-)/2 is the current velocity and u = {b^ — &-)/2 is the osmotic velocity. 
Further, we assume the velocities satisfy the Fokker-Planck equations 

-^ + V-(pv) = 0, u-DV\np = 0. (6) 

Now, we get equations that describe time dependence of u and v. 
du 



Z^ = -DV^v-V{v-u), 

^ = -Ivt;^ + ^^^Vn^ + (1 - 2si + 2s^)DW^u - -WV, 
at 2 2 m 



(7) 



here and hereafter, we assume the current velocity has no rotation. For a case Si = S2 = 0, 
they coincide with the Nelson's equations, for $2 = 0, they become the Davidson's equations 
[T2] . and for Si = s^ = 1/2, they become the equations we proposed in previous paper 
|13] . The second equation is derived by the sum of two equations ([1]), then there is another 
independent equation. The subtraction of Eqs. dl]) yields 

2(si - l){u ■ V)v + D{2si - 2s2 - l)V^v = (3v. (8) 

This equation will determine the magnitude of friction coefficient mp. When /3 = and 
V is constant, the equation becomes trivial. In general, right-hand side (RHS) may not 
be proportional to the current velocity, then, we shall need an extension /3vj — >■ l^ijVj by 
modifying Eqs. ([2]). 

3 Quantum process 

The original concept of Nelson's idea is to lead the Schrodinger equation from Eqs. ([7]). For 
convenience, we first consider the Lagrangian and the Hamiltonian from the equations of 
motion. According to [16], they are given by 

L = -mv^ + (1 - 2si)(Vlnp)2 + mD\l - 2si + 2s2)V^lnp- V, 

1 mD'^ ^^^ 

H = -mv^ {l-2si){V\npy -mD^{l-2si + 2s2)V^\np+V. 

When the Hamiltonian is constant written by E, the current velocity is zero, S2 = 0, and 
f{x) = ^J p{x), we get the equation 

E = -2{l-2si)mD'^—^ + V. (10) 

It looks like the Schrodinger equation if 2(1 —2si)mD'^ is constant and corresponds to fi^ /2m. 
If we use p instead of /, the first term of RHS corresponds to the quantum potential. Since 
p is real, the function / is also real. It approaches more to the Schrodinger equation if we 



consider the time dependence of p and v. Assuming v is gradient of some function, given by 
D'Vw with constant D', and S2 = 0, Eqs. ([7]) can be rewritten as 

^dOnpl ^ _jjjj'^2^^ _ DD'^^ . vinp, 

^^ (11) 

D'^ = -\v' + ^-^u' + (1 - 2s,)D'y\\np) - ^V. 
at 2 2 m 

Giving a new function / = g^^^^ ''+*"'", we have 

VV = (tV Inp ■ V Inp - a^Vw ■ Vw + -V^ Inp)/ + m(Vlnp ■ Vw + V^w)f (12) 

4r ^ 



If we take a = D' /2D^^1 — 2si, it becomes 



Then, this can be 



2(1 — ZsijU"^ ot m V ot 



2iy/l - 2simD^ = -2(1 - 2si)mDW + Vf, (14) 



Let D = h/2y/l — 2sim or si = 1/2 — h^/Sm'^D'^, it is the same to the Schrodinger equation 

Assuming /(t, a;) = f(x)e~'^^^^^ (i.e. aw(t, a;) = aw{x)—Et/h) for the case of constant energy, 
it is consistent with Eq. (fTOj) . A difference from quantum mechanics is that the function / is 
scalar while the wave function is vector in a Hilbert space. The extension to a Hilbert space is 
necessary to assign quantum charges and consider multi-particle system. This extension is not 
difficult if we employ infinite number of probability densities and velocities. Taking stochastic 
process x^^' for each component of a Hilbert space and giving mean velocities for each process 

(k\ (k) 

b, and 6_ , with k runs from 1 to infinity. For all the velocities, we apply the modified 



Langevin equations and the Fokker-Planck equations, i.e. Eqs. flMJ) and Eqs. Qj. Then we 
obtain the Schrodinger like equation for each k to the function /'^'^'^ = e^^^'^ '^*°'" . As for 
the potential V^^\ it is not necessary to take independent function in each component. That 
is, they can be functions of other components V^''^ = \/('^)(/(^), /(2), ■ ■ ■ ). When the potential 
can be written in linear combinations, we can take V^^'' fk — )■ "^Vkifi. More generally, the 
potential will be given by operator system with the function of fk through other components 
so that we may write V^''^ fk — >■ Vfk- Rewriting fk = ip, the equations turn into 

ih-^ = -— V^T/; + Vij. (16) 

ot 2m 

Now we can deal with any potential and any state of quantum mechanics as long as non- 
relativistic approximation is valid. Compared to the second quantization formalism, it is 
less convenient to calculate cross section, decay rate, branching ratio, and so on, due to 



the difficulty of formulating relativistic Brownian motion [T7] . However, when one considers 
a relationship between quantum and classical mechanics, above quantization scheme will be 
useful. As it is easy to contain classical equations, we can predict intermediate states between 
classical and quantum solutions to check the method. 

Note, as we focus on low energy physics, the system with a single component function 
is sufficient for the description of dynamics. For the moment, we do not mention about 
a Hilbert space. In addition, instead of dealing with many particle system with the wave 
function, we can employ macroscopic physics by reinterpreting p as the number density. In 
that case, Langevin equations should be modified but the same equations can be used with 
some simplification [18]. 

4 Diffusion process 

It can be inferred from Eq. ( IT4l) that the flip of the sign of 1 — 2s i reproduces the diffusion 
equation. Only, it makes a wrong sign compared to the actual equation. To produce proper 
one, we set /' = e'^^^P"'^ ^ and a' = D' l2D\J2s\ — 1 with assuming 2s\ — 1 > 0, then 

2^2si - ^mD-^ = 2(2si - l)mD'^V^f' - 2mD^S2{V^ In p)f' + Vf, (17) 

When S2 = and V^ = 0, it becomes 

df 



^ V2s^^lDV'f. 



Thus /' obeys the diffusion equation with the diffusion coefficient Dq = \/2si — ID. This 
function is not directly related to the probability density so that we need an additional 
equation. To compute the dynamics, we consider the function / = e2^^P~^°-^ which satisfies 

^ = -x/2i7^lDVV. (19) 

Giving initial conditions of / and /', the probability density p = f f can be computed. 
With a famous solution in two dimensions, i.e. / = foe^ /^^°*^*"*"*°''/47ri5o(^ + ^o) ^-nd /' = 
/oe"''^/^^o(*+*o)/47rL)o(t + t'^), the probability density becomes 

f — f' r^(tO-t'o) 

n = - f 4Do(t+to){t+to) (20) 

^ AnD,it + to)it + t',) • ^ ' 

In this way, time evolution differs from the diffusion equation in general. For a special case, 
when / is constant, the time-dependence becomes the same. The equations with the system 
of / and /' describe general dynamics of the diffusion equation. 

For the case / = const, or D'w = — ■\/2si — IZ^lnp, the probability density obeys the 
diffusion equation. This case means v is proportional to u, actually v = —y/2si — In. 
Suppose a solution p = e~^ ^^^"^/A-KDot, we can calculate the friction from Eq. ([8]). Since 
u = — a;/2\/2si — It, Eq. ([8]) reduces to /3 = — (si — l)/y/2si — It. This is one prediction for 



a particle obeying the diffusion equation in our system. It can be rewritten by observable 
parameters as 

.2 „,2 



/3 = -2DoW{-7rr^p)^----, (21) 



where W{x) is the Lambert W-function. For large t or small p, it can be approximated by 
P f» 2t[Dqp{u^ — v'^)/\u\\v\. The sign of the friction force depends on magnitudes of two 
velocities; it is positive when \u\ > \v\ and negative when \u\ < \v\. For general S2, the 
modified Nelson equations do not reduce to single equation in the above scheme. It becomes 
much harder to solve the equations. 

5 Soliton 

We have analyzed si > 1/2 and si < 1/2, and this section deals with si = 1/2. The 
case Si = 1/2 and S2 = illustrate just classical physics. This is natural because this 
limit corresponds to h ^ 0. Under the words of stochastic process, these assumptions lead 
l{d+ + d^)b± = ±l3v + K + dB. Then we get {d+ + d^)v = K + dB and \{d+ + d^)u = f3v. 
Since the average of {d+ + dJ) /dt corresponds to the Euler-Lagrange derivative, ordinary 
equations of motion realize. Precisely, they can be 

{- + v-V)v = K, {- + vV)u = l3v. (22) 

The first equation can lead Newton's equation of motion or it is related to the Navier-Stokes 
equations (to derive them we need to treat fluid density properly, see [16l [H]). The second 
equation is not trivial so it is a constrained system. Using the Fokker-Planck equations, it 
can be (n ■ V)v + /3d = 0. 

For the Hamiltonian with zero current- velocity, the condition si = 1/2 makes 

H = -2s2mD'^V^\np + V. (23) 

When the total energy and the potential are constant with written hy E — V = AE, it can 
be 

AE = -2s2mD^V^ In p. (24) 

In this section, we assume S2 > and AE > 0. For one spatial dimension, the probability 
density becomes Gaussian. For two dimensions with circular symmetry, the solution is p = 
(ci + C2r)e~^^^ /4.s2mD ^ y^iiQYe c\ and c^ are constant. For three dimensions with spherical 
symmetry, it is solved as 

p = c^(ri^^-^^l\ (25) 

When C2 > 0, it can be normalized by 



^r^ = \^ciGZ ( 



clAE 



48s2mD2 



|-1,0), (26) 



where G is the Meijer G-function. When the current velocity is non-zero constant, we have 



E mv 



V 






-V ■ Vlnp. 



(27) 



Writing E — ^mv"^ — V = AE, the solution can be written in the same way with p(t, x + vt): 



n ^ CtC 1232 "1-D2 ^/ V ^ ' ' 



{21 



This solution contains both of Gaussian soliton and bubble soliton. When C2 = it is a Gaus- 
sian and otherwise it can be seen as a bubble. The size of a bubble is tq = {6s2C2mD'^ / AEy^^ 
and the thickness of a film mainly depends on Ar = ^ Vls^vaD'^ j AE . The soliton will be 
stable as it is a general solution in that environment. By changing s\ and S2, the soliton will 
collapse in several ways. 

More generally, when v = VVw, f = e2 1'^/'+**"'^ and b = D' /\^8s^D^, we have 



^ M^ ^ ^^ 4s2D^ dt As2mD^ D' dt '^' 
where * denotes complex conjugate. Then we get 



iVS^mD^ = -4:S2mD^V'^f + {s2mD^{V\nf*fy + V)f. 



(29) 



(30) 



Thus the equations reduce to one differential equation for complex function / as in the same 
way of the Schrodinger equation. This equation will generalize the above soliton solution. 



6 Periodic localization 

The Schrodinger equation appears by taking S2 = and p = p in the Hamiltonian system. 
However, even when S2 7^ 0, we can get the same equation. By reparametrizing the function 
P = /2(i-2si+2s2)/(i-2si)^ ^Yie constant Hamiltonian without current velocity becomes 



E 



2(1 - 2si + 2s2)2mL)2 vV 



l-2si 
If 1 — 2si > and E — V > 0, a typical solution is 



/ 



/ = Co 



cos 



Axi 



Cl 



COS 



X2 

Ax, 



C2 



COS 



V. 



AX3 



For simplicity, we choose ci = C2 = C3 = and Axj = Ax, then 

Ax = 



'6(l-2si + 2s2)2mL)2 



[l-2si){E-V) 



In the view of probability density, it becomes 

l-2si 

r Xl T 2JJ 

P = Po 



COS 



Ax' 



2(l-2si+2s2) 



COS 



X2 



Ax- 



2(l-2s^+2s2) 



COS 



Ax- 



C3J 



(31) 



(32) 



(33) 



2(l-2s^+2s2) 



(34) 



equation (solution) 


Sl 


S2 


f 


quantum 


si = 1/2 - hySm^D^ 


S2 = 


complex 


diffusion 


si > 1/2 


S2 = 


real 


classical 


si = 1/2 


S2 = 


real 


soliton 


si = 1/2 


S2>0 


real 


periodic localization 


si < 1/2 


S2^0 


real 


NLS 


si = 1/2 - h'^/Sm'^D^ 


S2^0 


complex 



Table 1: Appearance of interesting equations and solutions are shown, depending on si, S2 
and /. If one believes that Planck constant is constant rigorously in quantum phenomena, 
Si should retain 1/2 — h'^/Sm'^D'^. However, if this is not true at some macroscopic scale, it 
can be any value within si > 1/2. 



If we take p as the number density, (1— 2si)/2(l — 2si + 2s2) ^ 1, and the Born-Oppenheimer 
approximation, it can roughly describe nuclei of solid state. Suppose the diameter of a nucleus 
is lOfm and the size of atom is 100pm, a realistic model can be obtained in the vicinity of 
(1 — 2si)/2(l — 2si + 2S2) ~ 10^. If Si is exactly equal to S2, strong localization is occurred 
by making the absolute values of Si and S2 very large. In other cases, it needs fine-tuning to 
keep 1 — 2si + 2s2 very small for the realization of such localization. 

Naively thinking, if it is regarded as a solid, phase transition can be performed in the 
framework of equations of motion. Varying the parameter si to Si > 1/2, it satisfies the 
diffusion equation so that nuclei in the solid will diffuse. Considering the diffusion process 
describe liquid state, this procedure can be interpreted as solid to liquid phase transition. 
Depending on how si changes, whole dynamics of this phase transition can be predicted. 

Next, let us consider generalization of above solution with non-zero current velocity. If 



we take si = 1/2 — nr /Sim?D^, S2 7^ 0, and / = ea 
Schrodinger equation: 



5 \np+iaw 



, a nonlinear term appears in the 



ih 



dl 
dt 



2m 



VV 



S2h'' 



2(l-2si)m 



{V'\nrf)f + Vf. 



(35) 



When the current velocity is zero and h is free, above periodic localized solution is appeared. 
The nonlinear term vanishes when the function is the plane wave, otherwise it affects to the 
energy at the same order of the kinetic term. By extending to a Hilbert space, the nonlinear 
term exists in each component, i.e. 



ih 



dfk 
dt 



h^ 

2m ■"" 



S2h' 



2(l-2si)m 



i\/Hiirjk)h + Vfk 



(36) 



where k is the index of the Hilbert space. Therefore, our nonlinear extension is quite different 
from ordinary interactions of quantum mechanics. 

Table 1 summarizes the equations and solutions that we have considered. In general, we 
need to solve nonlinear differential equations ([7]), but when 1 —2si > and v = D'Vw, they 
reduce to one complex equation. For the case 1 — 2si < and S2 = 0, the equations become 
two independent diffusion equations about / and /' to give the density p = ff. A remaining 



8 



region (1 — 2si < and S2 7^ 0) is much more difficult to be solved since the equations cannot 
reduce to single equation. 

7 Conclusion 

We have proposed a general framework that includes the Schrodinger equation, the diffusion 
equation, the soliton equation, and the Newton's equations of motion. Based on the modified 
Langevin equations for non-differentiable processes, the system can describe various motions 
depending on two parameters si and S2. Keeping S2 = 0, the parameter si can connect 
quantum, classical, and diffusion equations. On the other hand, taking non-zero value for S2, 
there exist some interesting solutions like soliton and periodic localization. 

Since the framework can provide a lot of predictions, the con&mation of the theory seems 
possible. For instance, it includes generalized diffusion equations with the system of p = //' 
and the new type of soliton equation with the potential V^ In p. Concretely, the potential 
terms u"^ and V ■ u with arbitrary coefficients are the prediction of the theory. If one finds 
such potential in soliton, diffusion, and quantum states, these states can be related in our 
simple model. 
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